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Abstract 

We link observational parameters such as the deceleration parameter, the jerk, the kerk (snap) 
and higher-order derivatives of the scale factor, called statefinders, to the conditions which allow 
to develop sudden future singularities of pressure with finite energy density. In this context, and 
within the framework of Friedmann cosmology, we also propose higher-order energy conditions 
■ which relate time derivatives of the energy density and pressure which may be useful in general 
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The current observational data from type la supernovae favours the cosmological 
models of dark energy of the phantom type [2]. Phantom is the matter which has a very 
strong negative pressure which violates the null energy condition g + p > (p - the pressure, 
g - the energy density) and, consequently, all the remaining energy conditions 3]. It has very 
severe physical consequences since it leads to classical and quantum instabilities "\. It is 
even argued that it should not be accepted and that there should be some other explanation 
of the phenomenon of dimming of supernovae ^]. On the other hand, the observational 
data is interpreted only in terms of the simple framework of Friedmann cosmology with an 
additional restriction of the barotropic form of the equation of state p = wg (w = const.) 
which links the pressure and the energy density of matter. However, in general, the equation 
of state does not necessarily have to be of barotropic type. The most interesting possibility 
is that it changes with time during the evolution of the universe which makes the barotropic 
index w = w(t) to be the function of time 



From the observational point of view one is able to expand the equation of state in series 
around the present moment of the evolution by taking appropriate derivatives of the pressure 
with respect to the energy density. This procedure, however, requires the usage of the time 
derivatives of the scale factor up to suitable orders which would show the "tendency" of 
the universe to change its equation of state in future or study its evolution in the past. 
Due to this, the observational parameters made of the time derivatives of the scale factor 
were named statefinders [8]. Of course, in a general case, the spatial derivatives of all 
the kinematic quantities in non-isotropic models may also be involved, but in this paper 
we restrict ourselves to isotropic models only. In fact, statefinders generalize such well- 
known observational characteristics of the expansion as the Hubble (first-order) and the 
deceleration (second-order) parameters. The third-order characteristics was called jerk ^] 
while the fourth-order characteristic was called either snap or kerk Q. In this paper we 
discuss even higher-order characteristics with the suggested names "lerk", "merk" etc. with 
the logical numbering them alphabetically with an addition of the term "erk" [lfl| . 

Phantom models of the universe (w < —1) admit a new type of singularity which is 
called Big- Rip (BR) |2j. At Big- Rip the energy density and pressure diverge as a result of 
having the infinite value of the scale factor a(t) at finite time. This is different from the 
ordinary Big-Crunch (BC) singularity which leads to a blow-up of the energy density and 
pressure for the scale factor approaching zero value at finite time. Due to time symmetry 
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the Big-Rip time singularity could have also appeared in the past but it is of course different 
from ordinary Big-Bang (BB) singularity. 

However, if one drops the assumption of an exact equation of state in the field equations 
one is able to face another t ype o f singularity in the universe, which is called a sudden future 
(SF) singularity [llj, |l2j, [l3|, |l4(. This is a singularity of pressure only, with finite energy 
density of matter. Unlike phantom, it appears for the matter fulfilling the strong energy 
condition g > 0, g + 3p > 0, the weak energy condition p > 0, g + p > 0, though violating the 
dominant energy condition g > 0, — g < p < g jl5, 16, 17 1 . The nature of SF singularities in 



terms of geodesic completeness has also been discussed 181 1 . 

nn n 

In this paper we discuss SF singularities [11J, 112J and generalized SF singularities dyj and 
relate their emergence to observational parameters composed of higher-order derivatives of 
the scale factor. These, on the other hand, may give some insight into the problem of the 
(changing in time) equation of state of the matter in the universe, since they are directly 
related to higher-order time derivatives of the energy density and pressure according to the 
expansion 

dp 1 d 2 p 

P = P ° + 1g '° ^ ~ 6o ^ + 2! dg 1 '° ^ ~ + °[( g ~ 

= Po + —{Q- Qo) + 7^ ~ (q-Qo) +0[{q-Qo)\ , (1) 

£o 2! g 6 

where index "0" refers to a quantity taken at the current moment of the evolution t = to- 

The SF singularities appear in the simple framework of Einstein-Friedmann cosmology 

governed by the standard field equations (we have assumed that 8nG = c = 1) 

fa 2 K\ 

( 'a a 2 K\ 

together with the energy-momentum conservation law 

Q = -^(Q + P) , (4) 

a(t) is the scale factor, K = 0, ±1 is the curvature index, g is the energy density, and p is 
the pressure. 

From (J3J) one can see that the singularity of pressure p — > oo occurs when acceleration 
a — > — oo, while a and a are regular. This can be achieved for the scale factor of the form 

a(t) = A+(a s -A)(f\ m - A (l - ±Y , (5) 



where a s = a(t s ) with t s being the SF singularity time and A,m,n = const. It is obvious 
from (jSJ) that a(0) = and so at zero of time a BB singularity develops. For the sake of 
further considerations it is useful to write down the higher-order derivatives of the scale 
factor i.e., 

a= -(a s -A)(iy" + A^(l-i) n " 1 , (6) 

m [m — 1) „ / t 

a s - A) 
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= 2 {as ~ A) \J. 

+ A n(n-l)(n-2) / _ tV" ? 
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or, for a general time derivative of an order r: 

(r \ m(m — l)(m — 2)...(m — r + l), / 1 

= * (°-- A )(r 



s , 

n— r 



( 1)r _ M n(n-l)(n-2)...(n-r + l) / _ (g) 



The main point is that the evolution of the Universe, as described by the scale factor ©, 
begins with the standard BB singularity at t — and terminates at SF singularity for t = t s 
provided we choose 

1 < 77 < 2, < 771 < 1 . (10) 

For these values of n and m the scale factor (J3J) vanishes and its derivatives ©-© diverge 
at t — leading to the divergence of g and p in 0-© (BB singularity). On the other hand, 
the scale factor (|SJ) and its first derivative (jBJ) remain constant while its second derivative (JZJ) 
diverge leading a d.vergeuee of pressure in Q only with finite energy density 0. It has 
been shown [l6| that more general sudden future singularities appear, provided we choose 

N < n < N+ 1, < m < 1 , (11) 

instead of (jlUj) . It means that for any integer N + 1 = r we have a singularity in the scale 
factor derivative (cf. (JUJ)), and consequently in the appropriate pressure derivative p( r ~ 2 \ 
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This, for any r > 3, gives a sudden future singularity which fulfills all the energy conditions 
including the dominant one 



Now, we introduce the higher-order characteristics or statefinders |8j, |3, |10( and relate 
them to the emergence of sudden future singularities. We also formulate the higher-order 
energy conditions which may be useful in general relativity. 

The well-known characteristics of the universe expansion are: the Hubble parameter 

H = l (12) 

and the deceleration parameter 

1 a da 

1 = —EH- = -72 ' ( 13 ) 



H 2 a a 2 

while the new characteristics are: the jerk parameter 



H 3 a d 3 ' 

and the "kerk" (snap) parameter [9|,[l 



1 a a a 2 

J = = -^sr • ( 14 ) 



H 4 a d 4 
the "lerk" parameter 

1 a^ a^a 4 



1 "a "a a 3 , , 

k = -771- = --T- , (15) 



(16) 



H 5 a d 5 ' 

and "merk", "nerk", "oerk", "perk" etc. parameters, of which a general term may be 
expressed as 

S W = (_!)H-i A a_ = ( _ 1)l+ iaa 

v ' H l a v ' a l+1 v y 

and its time derivative reads as 

= H[i(q + l)x^- (x {i+1) + x^)] . (18) 

A possible blow-up of statefinders may easily be linked to an emergence of singularities. 
In particular, this may be the signals for SF singularities. In fact, it is possible to formulate 
higher-order energy conditions which may be related to statefinders. We will study this in 
what follows. 
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The application of the definitions of the parameters (fT2*|l , (fT5j) , H3J) , (|ToT) , and (JTHJ) gives 
the following equalities for the time derivatives of the Hubble parameter [10|: 



H 
H 
H 
H 



-H 2 (q+1) , 

H 3 (j + 3q + 2) , 

-H 4 [k + 4j + 3q (q + 4) + 6] , 

H 5 [l + 5k + 10 j{q + 2) + 30q(q + 2) + 24] 



(19) 
(20) 
(21) 
(22) 



Using the equations (0), (j3J), (JU), (fT2"|) and (fTSj) the three energy conditions (weak, strong, 
and dominant) are equivalent to 

K 



g + p = 2H z (q + l + 



a 2 H 2 



g - p = 2H* \-q + 2 + 2 



ii 



a 2 // 2 



> 0, £ > , 



> 0, £ > 



e + 3j9 = 6g# 2 > 0, g > . 

In fact, the weak energy condition requires both 
application of (J2J), (J3J) and (J1J gives 



(23) 

(24) 
(25) 



and 



K 1 / 



to hold. Notice that the 



(26) 



a- 2 

which, after further differentiation with respect to time, and the usage of the formulas (|19j) - 
()22|) , may give the higher-order weak energy conditions which relate time derivatives of the 
energy density g and the pressure p. These are 

K 



g + p 



2H 6 [j + 3q + 2 + 2 



a 2 // 2 



> 



e > o 



g + p = 2H 4 [k + 4j + 3g(g + 4) + 6 



+2^^(4 + 3) 



a- 



!#2 



> , £ > 



'g + p = -2H 5 [l + 5k + 10j(q + 2) + 30q(q + 2) + 2A 
K 



+2 



a 2 H 2 



(j + 9q + 12) 



> 0, Q > 



(27) 



(28) 



(29) 
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The application of 



and (|27|) allows to write the time derivatives of pressure 

K 



P 

P 
P 



2E A 



2H 



-3 + 1 + 



a 2 H 2 



k + j-3(q+l)- 



K 



"(9 + 3) 



a?H 2 

-2H 5 [l + 2k + j(q - 1) - 6q 2 - 21q - 12 
K 



- (j + 9g + 12) 



l H 2 



(30) 
(31) 

(32) 



From (jHDl), (|3*T|l and (|3~2"|) one can see that the consecutive SF singularities as given by the 
condition (|TT|) are strictly related to a blow-up of the statefmders. In fact, for N = 1 there 
is blow-up of the deceleration parameter q (p diverges), for N = 2 there is a blow-up of the 
jerk j {p diverges) etc. 

Now we can write down the higher-order dominant and strong energy conditions, i.e, 

K 



Q-P 
Q + 3p 

Q-P 

Q + 3p 
"g - p 



-2H 6 [-j + 3q + 4 + 4 



a 2 H 2 



> 



Q > 



-QH 3 (q + j) > , g>0 
2H 4 [-k + 2j + 3q(q + 6) + 12 

a z H z 

6H 4 [k + 2j + q(q + 2)] > , g > 

-2H 5 [—l + k + 2j(iq + 11) + 6(7q 2 4 
K 



17 q + 8) 



+4 (j + 9q + 12) 



l R 2 



> , g > 



'£ + 3p = -QH 5 [I + 3k + 2j(5q + 3) + 6q(q + 1)] > 0, Q > 



(33) 
(34) 

(35) 
(36) 

(37) 
(38) 



Referring to SF singularities, one can see that it is not possible to fulfil any generalized 
dominant energy condition if any of statefmders q,j,k,l etc. is singular. This is because 
in the appropriate expressions (J2l?|) and (|24|). ()27|) and etc. the signs in front of corre- 
sponding statefmders are the opposite. 

Having given the higher-order time derivatives of the energy density and pressure related 
to statefmders, one is able to invent some more sophisticated energy conditions which may 



be obeyed by presumably singular cosmological models 
energy condition 
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3- 



As an example consider an 



ag > p , 



(39) 
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with a = const., which, after the application of (J2J) and (|3U|) gives 



i>|l-^fl + ^l. (40) 



2H J \ a?H 2 

and this may prevent the emergence of SF singularity for N = 2 in ([lljh 

Finally, the observational determination of the value of jerk by using type la supernovae 
sample has been performed ^| and it claims that jo > 0. Similar investigations may 
perhaps also be possible for higher-order characteristics such as kerk/snap, lerk etc. Some 
hints about that are given in M. 
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